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(Dated: February 4, 2008) 

Using 65 million Y(4S) -> BB events collected with the BABAR detector at the PEP-II e+e" 
storage ring at the Stanford Linear Accelerator Center, we measure the color-favored branching 
fractions B(B° -» D+tt~) = (2.55±0.05±0.16) x 10~ 3 , B(B° -> D* + tt~) = (2.79±0.08±0.17)xl0 -3 , 
B(B~ -» D°7r") = (4.90 ±0.07 ±0.22) x 10~ 3 and B(B~ -> D*°7r") = (5.52 ± 0.17 ± 0.42) x 10~ 3 , 
where the first error is statistical and the second is systematic. With these results and the current 
world average for the branching fraction for the color-suppressed decay B° — > D'*'°7r°, the cosines 
of the strong phase difference 8 between the 1 = 1/2 and 1 = 3/2 isospin amplitudes are determined 
to be cos 8 = 0.872±o.oo?1o.om for the B -> Dtt process and cos 8 = 0.924±g;^?±g-.g|| for the 
B — > D*tt process. Under the isospin symmetry, the results for cos 8 suggest the presence of final- 
state interactions in the Dtt system. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



The B — > Dtt and B — > D* tt processes provide very 
good opportunities to test the theories of hadronic B- 
meson decays due to their clean and dominant hadronic 
decay channels. With the development of heavy quark 
effective theory (HQET) [Hi] and soft collinear effective 
theory (SCET) [3j,|j], the theoretical description for these 
hadronic decays has improved considerably, and the fac- 
torization hypothesis in heavy quark hadronic decay has 
been put on a more solid basis. The three decay am- 
plitudes A for B — > Dtt can be expressed in terms of 
two isospin amplitudes, A 1 / 2 and A 3 / 2 , under the isospin 
symmetry of the strong interaction: 

A(B° -> D+7T-) = s/Tj3A 3/2 + VV3A 1/2 , (1) 
V2A(B° -> D\°) = yA~/3A 3/2 - y/2/%A 1/2 , (2) 

A(B- - D°7T-) = V3A 3/2 , (3) 

where isospin amplitudes A-^/ 2 and A 3 / 2 correspond to 
the transitions into Dtt final states with pure 7=1/2 
and 1 — 3/2 isospin eigenstates 0, An identical de- 
composition holds for B — > D*tt decays. The isospin am- 
plitudes are not necessarily the same in the B — > Dtt and 
B — ► Z3*7r systems. In the context of QCD factorization 
0, A x i 2 and A 3/2 for B ^ Dtt (similarly for B — > £>*7r) 
are related by 

— ip- = l + O(A50DM), (4) 

where mj is the 6-quark mass and Aqcb is the QCD scale. 
The deviation of the ratio Ai/ 2 / {-\f2A 3 / 2 ) from unity is 
a measure of the departure from the heavy-quark limit. 
The QCD factorization implies that the relative phase 5 
of A1/2 and A 3 / 2 is 0(t1qcd / Tit) ■ Final-state interac- 
tions (FSI) in the 1 — 3/2 and 7=1/2 channels can lead 
to a non-zero 8. A large value of 5 will substantially sup- 
press the destructive interference for the color-suppressed 
decay B° -> D^°tt°, thereby increasing the associated 
branching fraction. 



Recent experimental results on the color-suppressed 
decay B° -> L>M°7r @, 0, @| provide evidence for a 
sizable relative strong interaction phase between color- 
favored and color-suppressed B° — > D^tt decay ampli- 
tudes. It has been suggested Q that improved measure- 
ments of the color-favored hadronic two-body decay of 
the B meson will lead to a better understanding of these 
QCD effects. Further experimental results on the color- 
favored decay B — > Dtt suggest the presence of final-state 
interactions in the B — > Dtt process This paper 

presents new measurements of the branching fractions of 
B- -> L>M°7r- and B° -> £>M+7r _ (charge conjuga- 
tion is implied throughout this paper) and of the relative 
phase 5. 

This analysis uses (65.2±0.7) x 10 6 BB pairs collected 
at the T(AS) resonance with the BABAR detector at 
the PEP-II asymmetric-energy storage ring during the 
2001-2002 data taking period. Charged tracks are de- 
tected by a 5-layer silicon vertex tracker and a 40-layer 
drift chamber. Hadrons are identified by measuring the 
ionization energy loss dE/dx in the tracking system and 
the opening angle of the Cherenkov radiation in a ring- 
imaging detector. Photons are identified by an electro- 
magnetic calorimeter. These systems are mounted inside 
a 1.5-T solenoidal superconducting magnet. 

Kaon and pion candidates are selected from charged- 
particle tracks using dE/dx and the Cherenkov light sig- 
nature. Each charged track, except the track used as 
the soft pion to reconstruct D* + — > D°tt + , is required to 
have at least 12 hits in the drift chamber and a trans- 
verse momentum greater than 100 MeV/c. D° and D + 
candidates are reconstructed in the K~tt + and K~tt + tt + 
channels, respectively. In each case, D meson candidates 
are required to have a mass within 3a of the mean re- 
constructed mass value, where the mass resolution a is 
approximately 7 MeV/c 2 for D° and 6 MeV/c 2 for D+. 
A vertex fit is performed on D° (D + ) candidates with 
the mass constrained to the nominal value [jjj]- A D° 
candidate is combined with a low momentum tt + or tt° 
to form a D* + or D*° candidate, where the tt° candi- 
date is formed from two photon candidates and must 
have an invariant mass between 120 and 145 MeV/c 2 . 
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Combinations with an invariant mass difference Am = 
between 143 and 148 MeV/c 2 for D*+ and 
between 138 and 146 MeV/c 2 for D*°, corresponding to 
±3cr about the Am peak, are retained. Each B meson 
candidate is reconstructed using the selected D or D* 
candidate and an additional charged track that is not 
consistent with the kaon hypothesis. 

To reject jet-like continuum background events, the 
normalized second Fox- Wolfram moment i?2 13], com- 
puted with charged tracks and neutral clusters, is re- 
quired to be less than 0.5. We also require | cos6>t| to be 
less than 0.85, where 9t is the angle between the thrust 
axis of the B candidate and the thrust axis of the rest of 
the event in the e + e~ center-of-mass (CM) frame. 

B candidates are identified using the beam-energy- 
substituted mass m,Es = \J (y^s/2) 2 — p* 2 and energy 
difference AE = E* — yfa/2, where E* and p* are the 
energy and momentum of the reconstructed B candidate 
and y/s is the total energy in the e + e~ CM frame. B sig- 
nal candidates have toes ~iBfl, the B meson mass, and 
AE ~ 0, within their respective resolutions. The resolu- 
tion in AE, (tab, for various B modes ranges from 15.7 
to 18.1 MeV. We require that \AE-{AE)\ < 3d ae- For 
events with more than one B candidate, a \ 2 is defined 
with the D mass mo, Am and their resolutions as 

/ m D - (m D ) 



X 



+ 



Am - (At 



(5) 



and the candidate with the smallest \ 2 1S chosen. 

The event yield n for each mode of B — > D^*'it~ is ex- 
tracted by fitting the ibes distribution of the selected B 
candidates with an unbinned extended maximum likeli- 
hood fit. The tties distribution is fit to the sum of a signal 
component, modeled as a Gaussian, and a background 
shape. The background shape is parameterized as the 
sum of a Gaussian, representing the peaking background 
events that peak in jiies, and a phase space parameteri- 
zation function [l4T ] representing non-peaking combinato- 
rial background and continuum events. The parameters 
describing the background shape, including the relative 
normalization of the peaking component, are determined 
by fitting Monte Carlo (MC) simulated samples, with the 
signal events removed. The total signal and background 
event yields, as well as the shape parameters describing 
signal events, are free parameters in the fit. The fitted 
771es distributions for each of the B meson decay modes 
are presented in Fig. [TJ The peaking background yield 
n p b is about (2-4)% of the observed B signal yield, as 
shown in Table Q] 

For each studied B decay mode of B — > Z)Wf, the 
branching fraction is calculated as: 



B(B -► L> w 7r) = 



2fN BS eB(D(*))- 



(6) 



Here N BB - is the total number of BB pairs; e is the ef- 
ficiency determined from signal Monte Carlo events; / 



represents /_) or /oo, the charged or neutral B meson 

production ratios at the T(4S I ), which we assume to be 
/+_ = /oo = 0.5; and B(D^) is the branching fraction 
of D or D* decaying to its reconstructed final state 
The branching fractions we obtain are reported in Table 

in 

The final states D^tt selected by this analysis are, 
in general, accompanied by some small amount of final 
state radiation (FSR). We model final state radiation in 
our experiment with PHOTOS [l5j |. which predicts that 
6-7% of our selected events, varying slightly with decay 
mode, are accompanied by an average FSR energy of 
about 17 MeV. Approximately two-thirds of this energy 
is produced in the initial B decay, while the remainder is 
generated in the decay. 

We summarize systematic uncertainties on the mea- 
surements from various sources in Table [TTJ AN BB - is the 
uncertainty on the total number of BB pairs in data. The 
error on the efficicency, Ae, is due to signal Monte Carlo 
sample statistics. The uncertainty from combinatoric 
background is estimated as the difference in the B yields 
obtained when fixing and floating the non-peaking back- 
ground parameters in the jjjes fit. The uncertainty from 
peaking background is estimated as the B yield change 
by varying the peaking background parameters and the 
ratio of peaking background to non-peaking background 
within their errors in the m-Es fit- The uncertainties due 
to the differences in D^*> masses and AE between data 
and Monte Carlo samples are estimated by comparing the 
efficiencies using their resolutions and means from data 
and Monte Carlo samples in the event selection. The un- 
certainty due to D vertexing is estimated by comparing 
vertexing perfomance in data and Monte Carlo samples. 
The uncertainties in tracking, particle identification, and 
7r° reconstruction efficiencies are due to potential residual 
inaccuracies in the Monte Carlo simulation, after correct- 
ing for known differences. The dominant uncertainty is 
from the branching fractions B(D^) and the track- 
ing efficiency. 

With the branching fractions of the four color-favored 
decay modes B° — > D^ + w~ and B~ — > D^°ir~ , as well 
as the two color-suppressed modes B° — ► D^°tt°, one 
can calculate cos<5. Following Ref. [l6[ (equations have 
been modified to use the notation from Ref. Q), cos 5 for 
B — > Dir (similarly for B — > D*tt) can be expressed as 



cos 5 



3T{D+n-) + T{D°n-) - 6T(D°n°) 



6a/2|Ai 



/2-A3/2I 



1-4 



3/2 1 



-AD *-), 



(7) 



(8) 



\A 1/2 \ Z = T(D + n-) + rCDV) - -T(D%~). (9) 

Using the measured branching fractions in this analysis, 
the ratio of the B lifetimes t(B-)/t(B°) = 1.071 ±0.009 
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5.2 5.21 5.22 5.23 5.24 5.25 5.26 5.27 5.28 5%?Vc 

mEs (GeV/c 2 ) 




FIG. 1: Fit of m E s distributions for the B -> D W 7r candidates in data: (a) £° -> D+tt~, (b) B° -> D* + tt~, (c) B~ -> D°tt~, 
(d) B~ -> D*°7T~. The fit is shown as a solid line and is described in the text. The background component (including peaking 
background) is shown as a dashed line. 



TABLE I: Yield of signal (n) and peaking background (n p b), efficiency (e), and branching fraction (B) for each B —* D^*'-k 
decay mode. 



Mode n n ph e (%) B (xlO -3 ) 



B° 


-> D+tT 


3593 ± 63 


114 ± 14 


22.8 ± 0.2 


2.55 ±0.05 ±0.16 


B° 


-> D*+ti- 


1411 ± 39 


28 ± 6 


30.2 ± 0.2 


2.79 ±0.08 ±0.17 


B~ 


^D°7v- 


4606 ± 70 


89 ± 14 


37.9 ± 0.2 


4.90 ±0.07 ±0.22 


B~ 


-> D*°w~ 


1297 ± 39 


51 ± 8 


15.5 ± 0.1 


5.52 ±0.17 ±0.42 



and the branching fractions B(B° — ► Z?°7r°) = 
(0.291 ± 0.028) x 1(T 4 and B{B° -> D* tt°) = (0.27 ± 
0.05) x 10" 4 [Hi, we calculate cos 5 and \A 1/2 /(V2A 3/2 )\ 
for B — > Z)7r and — > D*tt decays. 

To estimate the systematic error on cos 8 for B — > Z?7r 
(and, similarly, B — > D*ir), we use a Monte Carlo tech- 
nique ficj]. We simulate 10 6 experiments, varying the 
measured branching fractions, the used color-suppressed 
decay branching fraction, and t(B~)/t(B°) about their 
central values according to Gaussian distributions where 
their errors are taken as the sigmas of the Gaussian dis- 
tributions, to calculate the cos<5. The correlation of the 
systematic errors between the two color-favored decay 
modes in the cost) calculation is taken into account. We 



assume the errors are uncorrelated between the color- 
favored and color-suppressed modes. The statistical er- 
ror on cos d is estimated in a similar fashion, with only 
the statistical errors on the branching fractions of color- 
favored modes are used in the procedure. The resulting 
normalized distribution of cos 5, ie., the estimated likeli- 
hood function of cos 6, is obtained. Figure [2] shows the 
likelihood function of cos 6 from the described experi- 
ments in which both the statistical and systematic errors 
are taken into account. 

We define ±ler confidence interval of cos<5 as the in- 
tegral of its likelihood function over the region around 
the nominal value of cos 5, which is calculated from the 
central values of the branching fractions, to 68.27% (half 
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TABLE II: Relative systematic errors in the branching fractions of B — > D^ir decays from different sources. 



Systematic error 


B° D+-K- 


B° -» D'+TT- 


B~ -> D°tt- 


B~ -> D*°ty- 




1.1% 


1.1% 


1.1% 


1.1% 


B(D (t) ) 


3.6% 


2.0% 


1.8% 


5.0% 


A/ 


1.6% 


1.6% 


1.6% 


1.6% 


As 


1.0% 


0.5% 


0.5% 


0.7% 


Non-peaking background shape 


2.8% 


0.5% 


1.9% 


1.3% 


Peaking background shape 


0.4% 


0.4% 


0.3% 


0.6% 


Data/MC difference of mo, Am 


0.2% 


1.3% 


0.4% 


2.9% 


Data/MC difference of AE 


0.5% 


0.2% 


0.6% 


0.7% 


D~ and D° vertexing 


0.2% 


0.1% 


0.1% 


0.1% 


Particle identification efficiency 


2.0% 


2.0% 


1.5% 


1.5% 


Tracking efficiency 


3.2% 


4.9% 


2.4% 


2.4% 


7T° reconstruction efficiency 








3.0% 


Total 


6.3% 


6.2% 


4.4% 


7.6% 



(a) 



0.6 0.7 0.8 0.9 1 1.1 



1.2 1.3-1.4 

cos8 




0.6 0.7 0.8 0.9 



1.1 1.2 1.3„1.4 

cos8 



FIG. 2: Likelihood function (arbitrary unit in vertical axis) 
of cos S obtained from the ensemble of 10 6 Monte Carlo ex- 
periments described in the text for process (a) B — » Dn and 
(b) B -* D*tt. The shaded area in the plots is 68.27% of the 
total area. 





0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.3 1 



FIG. 3: Likelihood function (arbitrary unit in vertical axis) 
of Ar = \A 1 /2/V2A 3/ / 2 \ obtained from the ensemble of 10 6 
Monte Carlo experiments described in the text for processes 
(a) B —> Dn and (b) B — > D*n. The shaded area in the plots 
is 68.27% of the total area. 



below and half above the nominal value) of the total area. 
The results are 



ros x _ o 72 +0.008+0.031 
CObO — u.o<z_ 007-0. 029 



(10) 



for the B — > Dir system and 



cos<5 - Q24+ aol9 + - 063 

cusu — u.»^<±_ ,017-0 .054 



(11) 



for the B — ► D* tt system, where the first error is statis- 
tical and the second is systematic. These results corre- 
spond to \6\ = 29.2°tof-tlfo and \6\ = 22.5°±lf ±l£ , 
for the B — * Dtt system and the B — * D*ir system, re- 
spectively. By comparing the likelihood function integral 
of cos S in region [0,1] with the full range integral, we ex- 
clude cos S > 1 at a probability of 99.9% for the B — > Dn 
system and 85.7% for the B — > D*ir system. 
Similarly, we obtain 



Al/2 



V2A 



3/2 



— «r;r;+0. 015+0. 042 

— U.UOO_ . 014-0. 042 



(12) 



and 



Ax/2 



V2A 



3/2 



= 0.624 



-0.027+0.065 
-0.026-0.063 



(13) 



for the B — > Dw and B — > D*ir system, respectively, 
where the first error is statistical and the second is sys- 
tematic. The likelihood function from the simulated ex- 
periments, with both statistical and systematic errors are 
taken into account, is shown in Fig. [3] 

In summary, we have measured the branching fractions 
for the color-favored B° -> DW+tt" and B~ -> D^°tt- 
decays. Using these measurements together with the cur- 
rent world averages for B(B° — > D°ir°) and B(B° — > 
D*°ir°), we extract the cosines of the relative strong 
phase 5 in the Dtt and £>*7r systems, and the ratios of 
the / = 3/2 and 1 = 1/2 isospin amplitudes. Our re- 
sults for the B — > D^tt branching fractions, except for 
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D*°tt , are consistent with the current world av- 



B~ 

erage values 12[ but have a better precision. The branch- 



ing fraction of B~ — > D*°tt~ from this measurement is 
greater than the world average by about 2a. Our results 
for cosS differ from unity by about 4.3<r for B — > Dir 
decays and 1.1a for B — > Z)*7r decays. The result of cos<5 
for B — ► Dir decays is consistent with the result in Refs. 
0, E3] , and under the isospin symmetry it suggests the 
presence of final-state interactions in B — ► Dir decays. 
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